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equation (including a new algorithm for reconnection processes) , and we present the results of calculation 
of some characteristics of a vortex tangle such as the total length, distribution of loops in the space of 



Abstract 

We present the results of simulation of the chaotic dynamics of quantized vortices in the bulk of 
superfiuid He II. Evolution of vortex lines is calculated on the base of the Biot-Savart law. The dissipative 
effects appeared from the interaction with the normal component, or/and from relaxation of the order 
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use the Langevin approach to the problem. In the present paper we require the correlator of the random 

force to satisfy the fluctuation-disspation relation, which implies that thermodynamic equilibrium should 
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■ be reached. In the paper we describe the numerical methods for integration of stochastic differential 
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CO , their length, and the energy spectrum. 
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t3 ■ 1 Introduction 
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, Quantized vortices appearing in quantum fluids influence many properties of the systems. In general, 

the set of these vortices represents a chaotic vortex tangle (VT) consisting of separated vortex loops 
(closed vortex lines). To describe the influence of VT, it is necessary to know the statistical description 
of VT at various moments of time evolution because at different processes quantized vortices reveal the 
different properties: thermodynamical equilibrium and nonequilibrium (turbulent). For example quasi- 
equilibrium features are essential in the process of fast quenching. In turn, in experiments with thermal 
flows and/or counterflows in He II a set of the vortices shows very nonequilibrium (turbulent type) 
properties. The statistical descriptions of these two cases are strongly different. Thus the distribution 
of vortex loops on their length n{l) — dN{l)/dl is governed by the formula [1] n(l) ~ l~ 5 ' 2 in the case 
of the thermal equilibrium, while in turbulent helium we have [2] n(l) ~ Z~ 4 / 3 . There exist many 
works devoted to numerical investigations of vortex tangles and to turbulent state of helium, e.g. [3- 
10]. Let's note that the mentioned calculations have been done in the local approximation framework. 
Originally simple vortex structures (VS) with time turn into a very strongly tangled system, if a self- 
crossing of the filaments happens in this system, the reconnection of vortex line occurs and thus the 
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vortex loops divide or confluence. Reconnections change the topology of vortex structures and affect 
the evolution of a VT. In the papers mentioned above, reconnection was simulated from the condition 
of the equality between the local and nonlocal contributions to the velocity of a vortex filament point, 
i.e., v n i w k/2ttA w vi w c/dn (R/ao)/4irR, where A = 2R/[c\n(R/ a®)] is the minimal distance between 
the pair of vortices, k is the quantum of circulation, R is the radius of curvature at the given point, c 
is the constant (~ 1), and a>o is the cutting parameter concerning with the radius of the vortex core. 
Thus, in those works only the distance between the points of a vortex line was chosen as the criterion 
for reconnection. In our opinion, this is slightly incorrect approach, because the elements of filaments 
can go away from each other and the crossing may not occur. In contrast to the mentioned above 
papers, we consider the entire Biot-Savart equation. Moreover, we take into account possible random 
disturbances in the system of vortices. The disturbances are simulated by addition of a new term to 
the Biot-Savart equation. The details see below. This statement is conventional for the description of 
dynamical systems with stochastical perturbations. Moreover, in this work the following condition of 
reconnections is proposed: if the elements of a VT have intersected during the temporal step of the 
calculation, the reconnection occurs. 



2 The problem statement and the dynamical equation 

We consider the dynamics of vortex loops in three-dimensional infinite space. The induced velocity of 
helium at a point r is defined by the Biot-Savart law: 

(S - r) x ds 



S - r| 3 



The formulae for the velocity of vortex line points, while without dissipation, takes the form [3]: 



^ = s 4/(w + fh ^ s , xS ,, (1) 



dt 4tt J |Si - S| 3 4tt I e V4 . ao 

where S(£, t) is the radius- vector of the vortex line points, £ is a parameter, in this case the arclength, S' 
is the arclength derivative, S + and 5_ are the lengths of two adjacent line elements that hold the point s 
between, and the prime denotes differentiation with respect to the arclength £. The second term of Eq.(l) 
is the local part of the velocity the first term is the nonlocal part obtained by integration on the rest of 
the vortex line and on all of the other loops. Taking into account the frictional force of vortices and the 
normal component of helium and the rapidly fluctuating random term (Langevin's force) , we obtain the 
equation for the dynamics of a point of the vortex line: 

S=S +a(s'x (v n -So)) -a'S'x [s' x (v„- So)] +A(t,£) (2) 

where (A{t,Q)=0, (^(ti,£i)A,-(t 2 ,&)) = D 5^5 (ti-t 2 )<f = D <5 y 5 (n(t) n{t')) , 

i,j are the spatial components; t\, ti are the arbitrary time moments; £i, £2 define any points on the vortex 
line; D is the intensity of the Langevin's force; a, a' are the friction coefficients; n(t) is the Gaussian 
white noise with (n(t) ) = 0, (n(t) n{t')) = 5 (t — t'). Let's assume further that the difference between 
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the normal and the superfluid velocities of helium equals to zero v„ = and neglect the term with a' . 
This statement correspond to absence of a heat flow. Thus, finally we obtain the dynamical equation for 
a vortex line: 

S=So -aS'x So +a(€)n(t), (3) 
In the integral form Eq.(3) is as follows: 

t t 
S(t, == S(to,0 + J B(£, t)dt + tr(0 J dW, (4) 

to to 

t 

where B(£,t) =So — aS'x So, and W{t) — J n(t')dt is the standard Wiener process. In our model the 

to 

initial conditions is six completely symmetrical rings with orientation making the total impulse of the 
system equal to zero. 



3 The numerical algorithm and the description of reconnections 

The Eq.(4) was solved by the Euler method: 

S n +i = S n + hB(£, t n ) + Vha(£)r] n , 

where S n is the approximate solution of the equation in a mesh point on time t n , h is the integration 
step on time in a mesh point t n , {r/„} is the consistency of independent between itselves normal random 
vectors with independent between itselves components in the aggregate r] n j (j = 1,2,3), having zero 
expectation value and the dispersion is 1. The component of vector r\ n was calculated by formula: 
rj n j = V~ 21nai cos(27ra 2 ), where ct\ and ct2 is random numbers from the interval (0,1) obtained by a 
pseudorandom-number generator. The Euler method is the first order on the mean-square approximation 
in the time step. The functions S',S",So were calculated as in paper [3]. To keep the calculation 
procedure coherent, points on the vortex line were added and removed was carry out as in the work [10]. 

The first step in the modeling of a reconnection process is selection of point pairs that are the candi- 
dates for reconnection. After the pairs was defined, it was assumed that the line segments between each 
of the pairs were moving with a constant velocity (according Vi,Vj) during the time step, as illustrated 
in Fig. 1. 

From the compatibility of equations 

Xi + V x ,ih + (x i+ i - Xi)sl = Xj + V x jh + {xj + i - xj)s 

Vi + V y ,ih + (y i+ i - y l )sl = y 3 + V y jh + (y j+1 - y 3 )s 
Zi + V z .ih + (z i+1 - Zi)sl = Zj + V z jh + (zj+i - z )s 
< si < 1;0 < s < 1 

it was determine the meeting of these line segments during the time step. Here (xi,yi, Zi, Xi + \, yt+i, z»+i); 
(xj,yj,zj, xj + i, yj+i, Zj+i) arc the coordinates of the first and the second pairs of points, accordingly; 
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Figure 1: The elements of vortex lines that must reconnect 



V Xt i,Vy t i,V Zt i;V x> j,V y> j,V z> j are the projections of the velocities of the points and the line segments on 
the coordinate axis. If the line segments had met, the reconnection occurs. Thus, if originally the points 
belong to the same loop, a pair of new loops was generated. Otherwise the confluence of the loops occurs. 

4 The results 

The initial radii of rings were R = 2 • 10~ 5 to . The initial condition was chosen in a way that the total 
impulse of the system was equal to zero. The rings were situated symmetrically at equal distance in pairs 
around the coordinate origin. The distance between them was d = 10~ 5 m. The parameters in Eq.2 are 
a = 0.0098, D = 4 • 10~ 5 m/s. 

Simulation was performed with a constant temporal step h = 5 • 10~ 8 s and initial steps along the 
vortex line was AZo = 2-7T • 10~ 7 m. The steps AZ along the vortex line was controlled later by the 
procedure of inserting and removing of points, so that AZo/2 < AZ < 2 AZo- It follows from Eq.4 that 
small (or those with a high curvature) loops move very rapidly and their dissipation (decrease in the 
size) is very high due to the friction. Therefore, small loops were removed during our calculations. Kinks 
appearing on vortex lines were removed also. For our case, the loops were cancelled if there were less 
than 5 points. 

Vortex configurations at various time moments are presented in Fig. 2. 

One can see a vortex structure with drastic evolution in time. After numerous time steps the system 
evolved into separated vortex tangles. It was noted that during the evolution, the vortex tangles arose 
and then vanished in different places. It resembles us the favor intermittency phenomenon in classical 
turbulence. 

From the simulation, several quantities were calculated and plotted as functions of time (see Fig. 3). 
One can see the phases of evolution for vortex structures. At the beginning, the total length, averaged 
curvature, and the vortex lines density increase at a steady volume. After a certain value of the density 
has been archived, many of small loops developers and the VT begin to decay. Later, one can observe a 
tendency to a fluctuating steady state. However, after the time t — 9.792 • 10 s, the volume suddenly 
begin to decrease. It is conform to disappearance of detached loops (see Fig. 2 c, d ). Then again the 
fluctuating steady state is reached. 
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Figure 2: The vortices at different times: a) t = 4.30 • 10~ 5 s, b) t = 9.792 • lCT 4 s, c) t = 1.0892 • 10~ 3 s, 
d) t = 1.1444 • 10~ 3 s. The configurations are plotted in a three-dimensional view. 



Figure 3: Plots of the total lines 1, the volume V, filed with vortices, and the vortex lines density L as 
functions of time. 
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We made also calculations for some statistical characteristics of a VS. The distribution of the loop 
number vs. its length and the VC energy spectrum were calculated. Figure 4 shows the distribution 
of vortex loops on their length at the time region t w 1.1443 • 10~ 3 s. Decreasing of the loops amount 
is described by the following function: n(l)dl ~ /-0.93 same dependance was observed both for 

different times and for each vortex tangle within the vortex structure. It is difficult to extract physically 
meaningful results from this data. On the one hand, the equilibrium has been reached. On the other 
hand, the distribution of vortex loops disagrees with the results obtained for thermal equilibrium and 
turbulent stationary state. 



Figure 4: The distribution of the vortex loops on their length at the time region t w 1.1443 • 10 3 s. 
The average kinetic energy of flow induced by a vortex loop can be evaluated as follows [11]: 

I L' L 1 00 

E= (J^ d3r ) = /^( 2w2 / /s(€i)S(6)e- tt(s « l) - s( « a) «id&\= jE(k)dk 

k * ' 

L V 

m = ^0rJi^JJ s (a) s fe) e " ik(s(6) " s(?2) ^^ 2 , 


where dflk = k 2 sm@kd@kd&k is the elementary volume, k is the wave vector, and p s is the superfluid 
density. For the isotropic case, the spectral density is expressed as follows: 

E[k) = W J J *|S(6)-S(6)| (dS(6) ' dS(6)) ' 


where k is the wave number. It is seen that there exist different regions of wave number k: small with 
respect to {V) 1 ^ , high with respect to 1/VL, and intermediate values. In the region of small wave 
numbers, E(k) oc k 2 , and for high numbers, E{k) oc kr 1 . 

Figure 5 shows our numerical results of the VC spectral density for the time t = 1.14445 • 10~ 3 .s . In 
the region of small wave numbers, simulation data fits the approximation E{k) oc k 2 . In the region of 
high wave numbers, it is difficult to extract a regularity from this data, since the spectral density and 
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the numerical error are of the same scale. In the intermediate region, the energy decreases according to 
the formula: E(k) oc fc~ 2 - 9 . 



Figure 5: The energy spectrum of the vortex tangle at t = 1.1445 • 10 3 s. 

5 CONCLUSION 

The obtained numerical results demonstrate that initially smooth vortex rings transform into a highly 
chaotic vortex tangle. In spite of that total length fluctuates about a constant value, we think that the 
thermal equilibrium state has not been reached yet. For instance, the spectral density of the energy 
agrees with the equipartition law only for a small k-zone. The vortex loop distribution on their length 
also differs from the one expected for a thermal equilibrium. This state is closer to the turbulent case 
observed by other authors. During evolution the vortex tangles arise and disappear in different places. 
It resembles the favor intcrmittency phenomenon in classical turbulence. Our preliminary simulation 
demonstrates that the Langevin approach applied in this paper is a very promising method for study of 
chaotic vortex structures. 

The work is supported by INTAS grant 2001-0618. 
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